Submitted to Phys. Lett. B 



E-Print hep-ph/980000 



Axion Emissivity from the Conversion of a Neutron Star into a 

Strange Star 

In-Saeng Suh * , Chul H. Lee t 

The Research Institute for Natural Sciences and Department of Physics, 
Hanyang University, 133-791, Seoul, Korea 
(February 1, 2008) 



Abstract 

The conversion of neutron matter into strange matter in a neutron star oc- 
curs through the non-leptonic weak- interaction process. We study the energy 
loss of the neutron star by the emission of axions in that process. Owing to 
that process, the neutron star will liberate the energy which can in no way 
be negligible as an axion burst. 
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I. INTRODUCTION 



The possibility of the existence of stable quark matter in the early universe or inside a 
neutron star or in relativistic heavy-ion collision experiments has been studied in many works 
[IJ. In this work, we consider the conversion of a neutron star into a strange star and its 
energy loss during the conversion process. In order to form strange matter in the interior of 
a neutron star, neutron matter should be converted into strange matter. Several conversion 
mechanisms have been discussed by Alcock et al. [f2|. When a neutron and a stable strangelet 
(strange quark matter droplet) meet, the neutron is readily absorbed, while a proton can 
coexist with a strangelet due to the Coulomb barrier between them ||. Therefore, a stable 
strangelet as a seed for the process will trigger a conversion in such a way that it grows 
by absorbing neutrons and, eventually, convert most of the neutron star into a strange star 
H|4[]. Then the conversion would liberate about 10 52 ergs in binding energy ||. 

A variety of mechanisms have been suggested for seeding the interior of a neutron star 
with stable strangelets so far @,||. They can be divided into two main categories ||; (1) the 
primary mechanisms in which the seed is formed inside the neutron star, (2) the secondary 
mechanisms in which the seed comes from the interstellar medium. In a viewpoint of the 
primary mechanisms ||, the central high densities and pressures in a neutron star may be 
sufficient for a phase transition from neutron matter to two-flavor quark matter to occur. 
Subsequently two-flavor quark matter can then easily decay into the lower energy strange 
matter through weak interactions. 

If we assume that strangelet has formed inside a neutron star, subsequent conversion 
process of the rest of the star is described as follows |3||| . The volume over which strange 
matter equilibrates was shown to be much smaller than that of the total strange matter 
region ||[7]], so that the problem can be treated hydrodynamically in one-dimensional ge- 
ometry. As the conversion front sweeps into neutron matter, the small region behind the 
conversion front has an excess of down quarks relative to strange quarks due to the flux of 
neutrons (udd) at the conversion front. The excess down quarks will convert into strange 
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quarks via the non-leptonic weak process ||, d + u — > s + u, as long as Hd > /^(^d > n s ). 
By this process a d quark can change itself to an s quark until the Fermi energies of all 
the flavors become equal and the energy per baryon drops in comparison to the ordinary 
two-component nuclear matter. Other leptonic decay processes are suppressed considerably 
||. The conversion will liberate ~ 10 52 ergs of energy (assuming that ~ 10 MeV is liberated 
per neutron converted). This energy will be radiated as neutrinos, photons, e + e~ pairs, and 
axions, etc. In this stage neutrinos and axions can escape the star. Many authors fUH have 
considered the escape of neutrinos. In this work we calculate the energy loss due to the 
emission of axions in the process of such a conversion of a neutron star into a strange star. 

At the chemical equilibrium between the quarks and the electrons, non- or semi-leptonic 
weak interaction is not important in neutrino and axion emission in quark matter because the 
weak interactions coupling s and u quarks are Cabibbo suppressed relative to the interactions 
coupling d and u quarks. For the nonequilibrium processes such as the conversion from two- 
flavor to three-flavor quark matter, however, as the strange quark semi-leptonic processes in 
neutrino emissions and non-leptonic processes in axion emissions are of significance. 

The axion is a pseudo-Goldstone boson which was introduced by Peccei and Quinn (PQ) 
to solve the strong CP-problem in a natural way [|TI]] . However, theoretical and experimental 
investigations give little guidance on the PQ symmetry-breaking scale, f a , and therefore on 



the mass of the axion [11|. The axion decay constant f a is related to the axion mass [12[] 

/0.62 x 10 7 GeV\ Jr 

m - K I To, J eK 

There are two generic types of invisible axions; the Dine-Fischler-Srednicki-Zhitnitskii 
(DFSZ) axion which couples to both quarks and leptons at tree level [[HJ and the hadronic, 
or Kim-Shifman-Vainshtein-Zakharov (KSVZ), axion which has no tree- level coupling to lep- 
tons but does couple to them at loop levels |T3|]. In invisible axion models, the axion mass 
m a is in principle arbitrary, however astrophysical and cosmological considerations |12j can 
provide an upper and lower bounds for m a . The astrophysical bounds on m a are due to the 
fact that axion emission is an additional energy loss mechanism for stars. 



If such axions could be copiously produced during the conversion of a neutron star into 
a strange star, it might drastically alter the energy budget of stars. The axions might thus 
deplete the stellar energy and change the usual course of stellar evolution. The emission 
of axions would have hastened the cooling process. The quarks are expected to be highly 
relativistic and degenerate in such matter with their Fermi energies in the range of 300 to 
500 MeV. A calculation of the axion energy flux from strange quark matter has been made 
by Anand et al. [ 15| . They have shown that the axion emission rate is several orders of 



magnitude smaller compare to the neutrino emission rate. 

We investigate here the energy loss of the neutron star during the conversion of non- 
strange quark matter into strange quark matter in the interior of the star. We assume 
that the axions escape freely out of the star as soon as they are produced, and that the 
temperature in the interior is constant during the conversion for the reasons that strong and 
electromagnetic interaction timescales are much smaller than those of weak interactions. 
Unless otherwise noted all equations assume h = c = k B = 1. 

II. AXION EMISSION RATE 

The axion is the hypothetical pseudo-Goldstone boson associated with the spontaneous 
breakdown of the PQ symmetry proposed to solve the strong CP-problem of QCD. We 
consider the invisible axions proposed by DFSZ that couple to both quarks and leptons at 
tree level. The interaction of axions with quarks can be written as [|12[ 

£int = {^j^j QWsQid^a, (!) 

where a is the pseudoscalar axion field, /„ is the axion decay constant. The CVs are model- 
dependent coefficients typically of order one. At tree level, while Cj = in the KSVZ axion 



model, they are nonzero for the DFSZ axion model. For example [16], 



C u = -2.76 - cos2/?, 

C d = -1.13 - cos2(3, (2) 



C s = 0.89 + cos2/3, 



where (3 depends on the Higgs sector. 

Since the conversion of a neutron star into a strange star occurs via the strangeness- 
changing non-leptonic weak reaction ||, the dominant axion emission process during the 
conversion period is 



u\ + d — > U2 + s + a . 



(3) 



For this axion emission process, there are four diagrams to be evaluated and one of them is 
shown in Figure 1. The matrix element from the diagram of Fig. 1 can be written as 



Sfi = u 2 



u d D F u s 
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y/2' 2 



7* (l- 7 5 ) s ine c 



ui, (4) 



where the subscripts 1 and 2 in u denote the initial and final u quark respectively, g is 
dimensionless weak coupling, q is the axion momentum and 8 C is the Cabbibo angle. The 
propagators are given as 

1 



iSp 



D, 



p 2 + fl-m q 

where p is the momentum of the quark and k is the quark four-momentum transfer p 2 — Pd 
carried by the weak boson. Since q pp, where p F is the quark Fermi momentum, and 
with the quark mass neglected, the matrix element, Eq. (4), simplifies to 



(5) 



Here we used the relation, for » k 2 [|T7 
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G f 



where Gf is the Fermi weak coupling constant and Mw is the H^-boson mass. For all 
processes the matrix element squared and summed over the color and spins is 



Sfi\ 2 T = 2 3 (jr) stn%cos 2 e c {2C 2 u + C 2 + C 2 S }\M[ 



(6) 



where 

\M\ = m 2 7m(1 - l5)ud u s j"(l - 7 5 )«i 
= Ps ■ P2Pd ■ Pi ■ 

The axion emissivity, the energy loss rate per unit volume and unit time due to the axion 
emission process, can be written as a multi-dimensional integral, 

£a = X 6)cfn i £ a |%|r 54 (Pl + Pd - P2 - Ps -Pa)S, (7) 

where i = Ui,d,U2,s,a and dUi = is the Lorentz-invariant phase-space element. 

The factor 6 stems from 2 spins and 3 colors, and | is explicitly shown in Eq. (7) to take 
into account that only left-handed helicity states of the u\ quark couple to W boson which 
mediates the transformation. Henceforth we neglect the masses of u, d, s quarks [§]. The 
statistical factor is 

5 = /i/ d (l-/ 2 )(l-/.). (8) 

Here the functions fi are the usual Fermi-Dirac distribution functions of the particle species 

h 

f= 1 

e t +1 

where //j are the chemical potentials usually of order of a few hundred MeV in the systems 
of interest. The distribution function depends on both the particle momentum and the 
species. The combination of distribution functions in Eq. (8) expresses the physical fact 
that for the process to occur, the initial u and d quark states must be occupied and the 
final u and s quark states unoccupied |)TB| . In quark matter, strong and electromagnetic 
interaction processes establish thermal equilibrium at a rate much faster than that of weak 
interaction processes. Therefore, we take the equilibrium Fermi-Dirac distribution function. 

In order to evaluate the energy loss rate Eq. (7), we first note that \q\ is small compared 
to the p^s and is therefore ignored in momentum conservation. Since highly degenerate 
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relativistic quark matter enables us to replace the momenta by their values at the top of 
the Fermi sea, the quark momenta |pj|'s are set equal to their respective Fermi momenta 
We shell neglect the effects of interactions on the particle density of states. Then, the 
phase space element d 3 pi = \pi\ 2 d\pi\dQi, where dfli is an element of solid angle, becomes 
d 3 Pi = \pi\ 2 dEidtti. 



In an integral over energies [|19(], we obtain 

Q Q 2 + (2n) 2 T 2 



J IlidEi5(Ei + E d -E 2 ~E S - E a )S = | 



er — 1 



where 

Q = E a + 5fi, 5/j = fj, d - /i s . 

We now integrate the above quantity over axion energies; 

2 + (2vr) 2 T 2 



T(T,5fi) = J dE a E 2 a Q- 



CT — 1 



«7--^ W <6) + A<S)] + ^rV. 



For an integral over angles [20 



A = j UjdQ^ipl + p d - p 2 - p s )\M\ 2 
1 (4vr) 5 



(9) 



P 3 f 2 

where we ignored the axion momentum, and j = Ui,d,U2,s. Therefore, Eq. (7) is written 

as 

£a = (^) {2C 2 U + C 2 + C 2 s }sin 2 e c cos 2 e c p* F A T(T, SfA). (10) 

Finally, from Eqs. (9) and (10), we obtain the axion emissivity 

£ a = 1.5 x 10 32 {2C 2 + C 2 + C 2 s }sin 2 e c cos 2 e c (^j ^-^^)Vj, eV (T, 8p) 

ergs/ cm 3 / sec . (11) 

where 



945 VMeVV 3tt 6 \MeW V Mel/ , 

11 / T \ 4 / ^ \ 2 
+ 90^ \MeVj \MeVJ ' ^ ' 

Figure 2 shows T MeV (Sfi) as a function of <5/i for the fixed temperatures T = 5, 10 and 
20MeV. For the chemical equilibrium, 8/i = 0, then T MeV ~ 6 x 10 4 for T = lOMeV. 
However, for the typical value of 5fi = 300Me1/, T Me v ^ 10 7 for T = lOMeV. 

III. RESULTS AND DISCUSSIONS 

Our main result for the total energy flux rate for the emission of axions from conversion 
of non-strange quark matter into strange quark matter is given by Eq. (12). Here we can 



carefully quote our bound on the axion coupling f a , using the simple inequality [21 



£ a ■ V ■ St < E m (13) 

The allowed energy loss, E m /5t is taken to be 10 52 erg/sec for definiteness, and the volume 
of axion emitting region V = f vri? 3 with R = 10 ~ 20£;m. As for the conversion timescale, 
if T = lOMeV, it varies roughly between one and ten minutes 0. So, we take the timescale 
typically St = lOOsec. Then, from Eq. (14), we can obtain the bound on f a , 

~ 10 5 GeV < f a . (14) 

We note that the bound based on Eq. (14) is legitimate only if the emitted axions never 
interact on their way out of the star. If the axion coupling is large enough, axions once 
produced may interact many times, namely be absorbed and reemitted more than once in 
the hot core of the star, and may be trapped thermally. 

As applications, Mikheev and Vassilevskaya [^] recently investigated the radiative decay 
of the axion a — > 77 in an external electromagnetic field in the DFSZ model. They concluded 
that the decay probability is strongly catalyzed by the external field, namely, the field 
removes the main suppression caused by the smallness of the axion mass. Therefore, the 



radiative decay of the axion in strong magnetic fields of the neutron star could give interesting 
astrophysical phenomena. Furthermore, we can consider that the sudden conversion from a 
neutron star to a strange star may account for the gamma ray bursts at comological distances 
H]. Most of the energy is probably released in the form of neutrinos. If a part of the total 
energy goes into 7 rays decayed from the axions, it will be large enough to account for the 
gamma ray bursts at cosmological distances and to explain their isotropic distribution. The 
outcome of such a conversion event will be the emission of as much as ~ 10 MeV of energy 
in ~ sec to ~ min as radiation with a typical temperature of tens of MeV. This conversion 
can be observed as a gamma ray burst j3|. 

In summary, as a possible mechanism for production of axions, we considered the conver- 
sion of non-strange quark matter into strange quark matter. We also estimated the energy 
loss of neutron stars through the emission of axions in addition to the cooling provided 
by the neutrino emission. During the conversion period, the important process is the non- 
leptonic weak-interaction. We assume that the axions are not trapped and escape freely 
out of the star as soon as they are produced, and that the temperature in the interior is 
constant during the conversion. We find that the energy carried away by axions during the 
conversion is not negligible. It is also found that the axion emission rate is three orders of 
magnitude larger compared to the chemical equilibrium case. In addition, we discussed the 
bound on the axion coupling. 
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Figure Captions 

Fig. 1 Feynman diagram for axion emission in the non-leptonic weak interaction. 

Fig. 2 Tuevi^lA as a function of b\i for the fixed temperature T = 10 MeV. The solid, dash, 
and dot lines correspond to T = 5, 10, and 20 MeV, respectively. 
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